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Abstract. Small three-dimensional strongly coupled clusters of charged particles in a 
spherical confinement potential arrange themselves in nested concentric shells. If the particles 
are immersed into a background plasma the interaction is screened. The cluster shell 
configuration is known to be sensitive to the screening strength. With increased screening 
an increased population of the inner shell(s) is observed. Here, we present a detailed analysis 
of the ground state shell configurations and configuration changes in a wide range of screening 
parameters for clusters with particle numbers N in the range of 11 to 60. We report three types 
of anomalous behaviors which are observed upon increase of screening, at fixed N or for an 
increase of N at fixed screening. The results are obtained by means of extensive first principle 
molecular dynamics simulations. 
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1. Introduction 

Coulomb crystal formation is among the most exciting cooperative phenomena in charged 
particle systems and has been observed in a variety of fields, including ultracold ions in Paul 
and Penning traps JHEIEI, electrons and excitons in semiconductor quantum dots Q and 
bilayers Coulomb crystallization occurs also in classical and quantum two-component 

systems such as electron-ion or electron-hole plasmas [171[8l or laser cooled expanding plasmas 
||9l , for a recent overview see [fTOll . Of particular recent interest has been crystallization of 
charged microspheres in complex plasmas in two dimensions ifTTl IT2l . as well as in three 
dimensions [13] since here the structure and dynamics of the individual particles is directly 
visible or recordable by standard CCD cameras, e.g. [,14] . 

From the theoretical side, the shell structure of spherically confined Coulomb crystals 
has been analyzed in great detail by computer simulations, e.g. [|T5l [T6ll and references 
therein. More accurate data including metastable states have recently been presented 
lfT7l[T8l[T9ll20ll2T| resulting in a very good understanding of these systems. However, in dusty 
plasmas the interaction of the two particles forming a crystal is screened by the surrounding 
electrons and ions which has a significant influence on the crystal structure. In [|22l l23l it 
was shown, by comparison with experiments, that the pair interaction is well described by 
an isotropic Yukawa potential and shell configurations for various values of the screening 
strength have been presented. It was found that screening leads to a cluster compression, a 
change of the average density profile |!24l [25l| and to an enhanced population of the inner 
shells. Yet a detailed understanding of how these shell occupation changes proceed is still 
missing. This is the goal of the present paper. 

Here we present a detailed analysis of the ground state configurations of mesoscopic clusters 
interacting via an isotropic Yukawa potential containing 11 to 60 particles in a wide range of 
screening parameters 0.0 < k < 20.0. While the general trend that with increased k particles 
move inward is confirmed, we observe several anomalies which are due to symmetry effects: 

1. upon K increase two particles move to the inner shell at once. 2. when the particle number 
is increased by one at a fixed k one particle move from the inner to the outer shell and 3. at 
very large k there exist cases of reentrent shell fillings: one particle returns from the inner to 
the outer shell. 

2. Model and Simulation technique 

We consider identical Yukawa interacting classical particles with mass m and charge q in 
a three-dimensionsal isotropic harmonic confinement potential described by the hamiltonian 

A'" N N 2 

H(r. V.) = i: + i: |r? + i: ^ . e—.. (1) 

This model has been found close to the experimental situation under which spherical dust 
crystals form ||26l . In the simulations we use dimensionless length and energy variables by 
introducing the units tq = {q^ /IQ-neaY^^ and Eq = (ag'^/327r^e^)^/^, respectively. 
This model was already used to find the ground state configurations and their energies for 
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Coulomb interaction in IfTTlfTSl . Here we extended the investigation to the ground states of 
finite Yukawa systems. To obtain the ground states we perform extensive molecular dynamics 
simulation using a standard simulated annealing technique, e.g. fi27ll . This is done by slowing 
down the particles by some friction in every time step, starting from a random configuration. 
A stable state is reached when the dimensionless force on each particle is zero (less than 10"^ 
in the calculations). It was observed previously for Coulomb systems that there exist different 
states with the same shell configuration, which differ with respect to the particle arrangement 
within the shells (fine structure) ifTTl [T9ll20ll . Here, these energy differences which are less 
than 10~^ in dimensionless units will not be resolved since this would blow up the whole 
analysis and we record only the energetically lowest shell configurations for a given value of 

Metastable states with a different shell configuration are sometimes energetically very 
close to the ground state. Also, their number is increasing approximately exponentially with 
N [[20| which requires special care in the numerical approach, in particular in the choice of the 
cooling speed. Also, for given parameters, the cooling process has to be repeated sufficiently 
often. In the present calculations, we typically used 10^- 10"^ independent runs for every set 
of (A^, k). While this does not guarantee that the true ground state is found it does ensure a 
sufficiently high probability that no other state with lower energy exists. As an independent 
tool to verify the results we performed for a number of cases standard Metropolis Monte Carlo 
simulations. 

We simulated particle numbers from 11 to 60 and screenings from k = to k = 5.0. The 
screening parameter was changed in steps of An = 0.1. When for some a configuration 
change at some critical k was detected, the calculation around this point was repeated with 
a substantial smaller k step to ensure an accuracy of ±0.05. The choice of this interval 
of screening parameters is motivated by the situation of typical dusty plasma experiments 
where k is around 1. Besides, it is of theoretical interest what will be the asymptotic shell 
configuration in the limit of a very short range interaction. To this end we also analyzed 
the ground state at k = 20 and recorded structural transitions occuring between k = 5 and 
K = 20. 

3. Numerical results 

3.1. Total energy 

A typical simulation result is shown in Fig. [T] where we plot the total energy per particle for 
the cluster A^ = 29 in the range of 0.0 < k < 5.0. As one can see the energy decreases rapidly 
with K by approximately one order of magnitude, due to the reduction of the pair interaction 
strength. The same behavior is observed for other particle numbers, as shown for A^ = 31 
in Fig. [2] and, for A^ = 57, in Fig. |3l Due to the exponential dependence on the distance one 
may wonder if the energy decrease with n follows an exponential law as well as is the case in 
macroscopic one-component Yukawa plasmas, e.g. [|28ll29l . 
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Figure 1. (colour online) Energy per particle of a Yukawa cluster with = 29 particles 
for screenings 0.0 < k < 5.0. The red solid [blue dashed] hne indicates the configuration 
(25, 4), [the configuration (24, 5)]. The vertical blue dashed line denotes the screening from 
which the configuration (24, 5) begins to occur in the simulations. The configuration (25, 4) is 
present in the complete range of screening. The inset shows the energy difference per particle 
of these two configurations in a small range of screening parameters around the critical value, 
where the ground state shell configuration changes from (25, 4) to (24, 5). The critical value 
is indicated by the vertical black dashed line in the figure as well as in the inset. 



The simplest fit for tlie ground state total energy per particle has the form 

and uses three k independent free parameters which are functions of the particle number. In 
the analyzed range of N this dependence is found to be close to A^^/^, for the two energies 
Eq and El, whereas the effective length ri in the exponent scales approximately as N^^^. 
Using the exact results for the ground state energies per particle from the molecular dynamics 



simulations we obtain the following best fit for the three coefficients: 

Eo{N) = 0.015 + 0.12iV2/^ (3) 

Ei{N) = - 0.81 + 0.92A^2/^ (4) 

ri{N) = 0.51 + 0. igiV^/l (5) 
In the Coulomb limit this fit reduces to 

eL{'^ = 0,N) ^ ^^^^^ ^ ^^^^^ ^ ^ iMN^/\ (6) 
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N 




Egs/N (MD) 


E^cs/N [Eq. m 


A(%) 


12 


0.0 


4.839 


4.656 


-3.8 


12 


4.0 


0.685 


0.736 


+7.4 


58 


0.0 


15.875 


14.788 


-6.8 


58 


4.0 


1.692 


1.902 


+ 12.4 



Table 1. Ground state Energies per particle from Eq. (O, compared to the exact results from 
MD simulations, and the relative error A, for some examples. 

This fit is useful to understand the main trends in the analyzed parameter range and reproduces 
the simulation data within several percent. Some representative examples are given in Tab.[Tl 
Further improvements can be easily achieved using e.g. the numerical results of ref. ||29l or 
the analytical expressions of ref. [|30il . but this is outside the goal of the present analysis. 

3.2. Structural transitions with screening 

The presented fit for ground state total energies Eq^ is a continuous functions of k and do not 
immediately reveal possible changes of the shell configuration. In fact, in many cases there co- 
exist several stationary states (shell configurations), the energies of which may become equal 
at a certain value of k. At this point a structural transition of the ground state is observed. 
This can be seen in Fig. [T]for the cluster with = 29 particles. For small k the configuration 
(25, 4) is the ground state until at the critical value of /t^ = 1.58 the configuration (24, 5) 
has the same energy and a smaller energy beyond this point, see inset of Fig. [U Thus, if k 
crosses k„. from below, one particle of the cluster moves from the outer to the inner shell. 
This ground state change is accompanied by a jump of the derivative of the exact ground state 
energy dEcs/dn at Kcr, so this structural transition resembles a first order phase transition. 

Figure |2] shows a more complicated example with two ground state changes occuring 
in a small range of screening parameters. For k < 1.5623 the ground state configuration is 
(27, 4) whereas at Kcri = 1.5623 the configuration (26, 5) becomes the ground state. Finally, 
at Kcr2 = 1.6142 this configuration is replaced by (25, 6) which remains the ground state for 
larger k. This behavior can be seen in the energy differences plotted in the inset of Fig. [2l 
Around the interval [ncri, Kcr2] all three states co-exist and have very close energies which 
illustrates the high accuracy and fine k— grid required in this analysis. 

These two examples are typical for most cases: at small k the cluster structure is strongly 
influenced by the spherical trap. In contrast, in the limit of very large screening the pair 
interaction tends to a hard sphere interaction and the clusters approach a closed packed 
structure. This is often a layered structure allowing for an optimal compression [|2T|. In 
between the two limits of long range and short range interaction the shell configurations 
change via one (or several) structural transitions where one particle from the outer shell moves 
to the inner shell as this configuration becomes energetically favorable. 

There are, however, several interesting exceptions to this general behavior. We observe 
three kinds of "anomalies" which will be analyzed in the following 
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Figure 2. (colour online) Energy per particle of a Yukawa cluster with 31 particles for 
screenings 0.0 < k < 5.0. The red solid line indicates the configuration (27,4) and the 
blue dashed [green dotted] line the configuration (26, 5) [(25, 6)]. The vertical dashed lines 
denote the beginning [blue for (26, 5) and green for (25, 6)] and the end [red for (27, 4)] of 
occurance of these configurations in the simulations. The inset shows the energy difference 
per particle for two stable states the red [blue] solid line for the configurations (27, 4) — (26, 5) 
[(26, 5) — (25, 6)] around the critical value of screening. The critical values for the changes 
in the ground state configurations are indicated by the vertical black dashed lines, both in the 
inset as well as in the figure. 

3.3. Anomalies of first kind: Correlated two-particle transitions 

Consider now the cluster N = 57, cf. Fig. [3l At small screening, the configuration (45, 12) 
is the ground state until at k^i = 0.10 one particle from the outer shell moves to the 
cluster center forming a new shell with the configuration (44, 12, 1). Thereby the second 
shell is not changed since it has a "closed shell" configuration with 12 particles. Besides 
this "normal" transition, at Kcr2 = 1-04 a new type of structural change is observed: the 
configuration changes according to (44, 12, 1) — > (42, 14, 1). This means, at this point a 
correlated intershell transition of two particles is observed. This unusual behavior will be 
called "anomaly of first kind". The reason of this anomaly is the particularly high stability 
of the closed shell configuration of shell two which dominates the structure up to rather large 
screening. In contrast, a configuration with 13 particles on the second shell is energetically 
very unfavorable, although it exists in a broad range of k values, in fact, the configuration 
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Figure 3. (colour online) Energy per particle of a Yukawa cluster with 57 particles for 
screenings 0.0 < k < 5.0. The red solid line indicates the configuration (45, 12) and 
the blue dashed [green dotted] line the configuration (44, 12, 1) [(42, 14, 1)]. The vertical 
dashed lines denote the beginning [green for (42, 14, 1))] and the end [red for (45, 12) and 
blue for (44, 12, 1), respectively] of occurance of these configurations in the simulations. The 
inset shows the energy difference per particle: the red [blue] solid line for the configurations 
(45, 12) - (44, 12, 1) [(44, 12, 1) - (42, 14, 1)] around the critical range of screening. The 
green solid line is the energy difference of the metastable configuration (43, 13, 1) to the 
current ground state, this configuration is never the ground state. The critical values for the 
changes in the ground state configurations are indicated by the vertical black dashed lines. The 
change (44, 12, 1) (42, 14, 1) at Kcr2 = 1-04 shows an anomaly of the first kind. 

(45, 13, 1) is never the ground state as can be seen in the inset of Fig. |3l 

The first occurence of an anomaly of the first kind is at = 30 where a transition 
(26, 4) (24, 6) is observed at k ^ 1.5. There is a total of 18 occurences of such anomalies: 
at iV = 30, 34, 36, 38, 40, 45 - 54, 57, 58, 60. The reason for this behaviour is that in all 
cases but for = 57, 58, 60 the new ground state configuration, e.g. (24, 6) at screening 
above n = 1.5, always forms a platonic body on the inner shell. This is a highly symmetric 
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Figure 4. (colour online) Ground states of small (11 < < 60) Yukawa balls for the range of 
screening parameter (0.0 < k < 5.0). The numbers on the bars denote the number of particles 
on the inner shell(s). The black circles indicate anomalies of the 1st kind. The white circles 
indicate the end of the screening range, where anomalies of the 2nd kind appear The ground 
states for a screening parameter k = 20.0 are plotted above the diagram for comparison in 
what range the ground states at k = 5.0 are stable. The cyan bar for TV = 43 at k = 20.0 
refers to a ground state of (11, 1) in the center region; it is the only time this configuration is 
part of a ground state. The dark blue circles just below k — 20.0 indicate anomalies of the 3rd 
kind, where a ground state configuration reappears with increased screening. 

configuration which obviously decreases the energy per particle better than by just adding 
one particle ifTTl [3T1l . For the cases = 30,34,36,38,40 the ground state configuration 
even change from one platonic body to another, while for the cases = 46 — 54 the system 
changes from 10 particles on the inner shell to the closed shell configuration with 12 particles. 
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The only configuration with 11 particles on an inner shell is found for the case = 43, 
at the very large screening value of k = 20.0, which leads to the conclusion that this 
configuration is energetically unfavorable. In the other three cases, = 57, 58, 60, the ground 
state configuration changes from (12, 1) to (14, 1) on the inner shells. Although a ground state 
configuration with (13, 1) particles in the cluster center is observed for some particle numbers 
in a certain range of screening parameters, the configurations (12, 1) and (14, 1) are far more 
often the ground state. 

These anomalies are shown in the full ground state diagram. Fig. IH by the black circles. 
The complete list is also shown in Table |2]by the bold numbers. 

3.4. Anomalies of the second kind: Reduction of inner shell population upon increase ofN 

Let us now consider changes of the total particle number at constant screening. The 
"normal" trend upon an increase of the particle number by one is, of course, that the new 
particle is added to one of the existing shells (leaving the other shells unchanged) or moves 
into the center opening a new shell. However, again, one observes exceptions from this rule, 
cf. Fig. m This effect was already observed for the Coulomb cluster {k = 0) with = 59 
[fTTl . It has the ground state configuration (46, 12, 1). Addition of another particle to the 
cluster gives rise to the configuration (48, 12). This is again a structural transition involving 
correlated behavior of two particles which we call "Anomaly of the second kind". In this 
particular case this transition is even associated with a change of the number of shells: the 
three-shell configuration [first appearing at A^ = 58] disappears again and, instead, a two- 
shell configuration is restored. This is, of course, a consequence of the particular stability of 
the latter which contains two closed shells with 12 and 48 particles, respectively. The closed 
shell configurations in dependence of the screening are given in Fig. [51 

While, in Coulomb systems, A^ = 59 is the only known case of an anomaly of the 
second kind, in Yukawa clusters this behavior appears quite frequently. The first occurence 
is at A^ = 11 for k values between about 2 and 4. Here, addition of a particle gives 
rise to the configuration change (10,1) — > (12,0), i.e. one particle moves away from 
the inner shell (the shell vanishes), and the population of the outer shell increases by 
two. There is a total of 20 such anomalous transitions observed for 18 particle numbers: 
A^ = 11, 20, 25, 26, 30, 33 — 42, 49, 55, 60. There are two particle numbers where this effect 
occurs two times: for A^ = 30 — > 31, in the k range [1.5623, 1.6142] the configuration 
changes from (24, 6) to (26, 5). Interestingly, for screening parameters just below this range, 
i.e. [1.4866, 1.5623] the inner shell loses even two particles, i.e. we observe the transition 
(24,6) — ^ (27,4). The second case where two such transitions occur is the transition 
39 — > 40. There for k between 0.2223 and 0.4179 the ground state changes according to 
(32, 7) — > (34, 6) whereas at k > 3.612 the configuration change is (31, 9) — > (32, 8). 

Finally, anomalies of the second kind which are additionally associated with vanishing 
of one "shell" (i.e. removal of one particle from the cluster center) are found 4 times: for 
A^ = 11 — > 12 the transition (10, 1) — > (12) is observed (see above). Return to a 
two shell configuration occurs three times: for A^ = 49 — > 50 we find the transition 
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Figure 5. (colour online) Shell closures for the ground state configurations for the first two 
shells in the range of screening 0.0 < k < 5.0. The particle number for the last closed 
shell is given by red dots [blue squares] for the shell 1 [2] in the considered range of screening. 
In some cases (k = [2.0 < k < 4.0] for shell 2 [shell 1] there exist different number of 
particles with closed shells, e.g. N = 57 and = 60 for k = 0. 



(36, 12, 1) — > (38, 12), for N = 55 — > 56, the transition (42, 12, 1) — > (44, 12) and, 
for = 59 — > 60, the transition (46, 12, 1) — ^ (48, 12)which is known from the Coulomb 
case (see above) and appears here in a narrow range oftmw /rvalues. The complete set of 
these anomalies is given in tabled 



3.5. Anomalies of the third kind: Reentrant shell transition upon increase of n 

Finally, there is a the third kind of anomalous behavior which deviates from the "normal" shell 
filling trend of increased populations of the inner shells upon increase of k at a constant A^. 
This tendency is never violated in the considered range of particle numbers, 11 < < 60, 
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N 


Her 


GS 


N 




GS 






GS 






GS 


11 


1.9 


(11,0) 


30 


1.5 


(26,4) 


42 


0.2 


(35,7) 


51 


0.5 


(41,10) 




> 5.0 


(10,1) 




> 5.0 


(24,6) 




1.1 


34, 8) 




1.1 


(39,12) 


12 


4.1 


(12,0) 


31 


1.56 


(27,4) 




1.8 


(33,9) 




> 5.0 


(38,12,1) 




> 5.0 


(11,1) 




1.61 


(26,5) 




> 5.0 


(32,10) 


52 


0.3 


(42,10) 


13 


> 5.0 


(12,1) 




> 5.0 


(25,6) 


43 


0.7 


(35,8) 




0.9 


(40, 12) 


14 


> 5.0 


(13,1) 


32 


0.6 


(28,4) 




3.6 


34, 9 




> 5.0 


(39,12,1) 


15 


> 5.0 


(14,1) 




1.5 


27,5) 




> 5.0 


(33,10) 


53 


0.2 


(43,10) 


16 


> 5.0 


(15,1) 




> 5.0 


(26,6) 


44 


0.7 


(36,8) 




0.7 


(41,12) 


17 


> 5.0 


(16,1) 


33 


0.2 


(29,4) 




1.5 


(35,9) 




2.9 


(40,12,1) 


18 


> 5.0 


(17,1) 




0.6 


(28,5) 




> 5.0 


(34, 10) 




> 5.0 


(39,13,1) 


19 


4.6 


(18,1) 




> 5.0 


(27,6) 


45 


0.5 


(37,8) 


54 


0.0 


(44,10) 




> 5.0 


(17,2) 


34 


0.2 


(30,4) 




1.2 


36,9 




0.6 


(42,12) 


20 


1.1 


(19,1) 




> 5.0 


(28,6) 




4.2 


(35,10) 




1.7 


(41,12,1) 




> 5.0 


(18,2) 


35 


0.3 


(30,5) 




> 5.0 


(32,12,1) 

\ ' 'J 




> 5.0 


(40,13,1) 


21 


1.2 


(20,1) 




3.2 


(29,6) 


46 


0.2 


(38,8) 


55 


0.5 


(43, 12) 




> 5.0 


(19,2) 




> 5.0 


(28, 7) 




0.9 


(37,9) 




1.5 


(42,12,1) 

\ ' 'J 


22 


0.3 


(21,1) 


36 


3.6 


(30,6) 




3.3 


(36,10) 




1.6 


(41,13,1) 

\ " 'J 




> 5.0 


(20,2) 




> 5.0 


(28,8) 




3.4 


(34, 12) 




> 5.0 


(40,14,1) 


23 


1.4 


(21,2) 


37 


1.7 


(31,6) 




> 5.0 


(33,12,1) 


56 


0.5 


(44, 12) 




> 5.0 


(20,3) 




4.2 


(30, 7) 


47 


0.5 


(38,9) 




1.0 


(43,12,1) 


24 


1.0 


(22,2) 




> 5.0 


(29,8) 




1.8 


(37,10) 




1.2 


(42,13,1) 




2.6 


(21,3) 


38 


2.0 


(32,6) 




2.1 


(35,12) 




> 5.0 


(41,14,1 




> 5.0 


(20,4) 




> 5.0 


(30,8) 




> 5.0 


(34,12,1) 


57 


0.1 


(45,12) 


25 


0.3 


(23,2) 


39 


0.2 


(33,6) 


48 


0.4 


(39,9) 




1.0 


(43,12,1) 




2.2 


(22,3) 




2.4 


32, 7) 




1.4 


(38,10) 




> 5.0 


(41,14,1) 




> 5.0 


(21,4) 




3.6 


(31,8) 




1.6 


(36,12) 


58 


1.0 


(45,12,1) 


26 


0.5 


(24,2) 




> 5.0 


(30,9) 




> 5.0 


(35,12,1) 




> 5.0 


(43, 14, 1) 




0.7 


(23,3) 


40 


0.4 


(34,6) 


49 


0.2 


(40,9) 


59 


0.5 


(46,12,1) 




> 5.0 


(22,4) 




> 5.0 


(32,8) 




1.1 


(39,10) 




0.9 


(45,13,1) 


27 


1.0 


(24,3) 


41 


0.1 


(35,6) 




1.4 


(37,12) 




> 5.0 


(44, 14, 1) 




> 5.0 


(23,4) 




0.6 


(34, 7) 




> 5.0 


(36,12,1) 


60 


0.1 


(48, 12) 


28 


0.1 


(25,3) 




0.8 


(33,8) 


50 


0.1 


(41,9) 




0.2 


(47,12,1) 




> 5.0 


(24,4) 




> 5.0 


(32,9) 




0.7 


(40,10) 




0.8 


(45,14,1) 


29 


1.6 
> 5.0 


(25,4) 
(24,5) 










1.6 

> 5.0 


(38,12) 
(37,12,1) 




2.9 
> 5.0 


(44, 15, 1) 
(43, 16, 1) 



Table 2. Table of structural transition points cf. Fig. |4]i. Bold values mark anomalies of 
the 1st kind, where the inner shell changes by 2 particles with increased screening {N fixed). 
The screening values displayed are the critical values (±0.05) up to which the configuration 
given in the 3rd column remains the ground state. 
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configuration 






11 12 


do 1) (12) 


1.9038 


4.0567 


20 21 


(18, 2) ^ (20, 1) 


1.0762 


1.1906 


25 — > 26 


(22,3) (24,2) 


0.2544 


0.5049 


26 27 


(22,4) ^ (24,3) 


0.7287 


1.0412 


30 — > 31 


(24,6) -> (27,4) 


1.4866 


1.5623 


30 31 


(24,6) (26,5) 


1.5623 


1.6142 


33 — > 34 


(28, 5) (30,4) 


0.2012 


0.2450 


34 — > 35 


(28,6) (30,5) 


0.2450 


0.3034 


35 36 


(28,7) ^ (30,6) 


3.1665 


3.6133 


36 — > 37 


(28,8) (30,7) 


3.6133 


4.1646 


37 38 


(30^7) (32,6) 


1.6679 


2.0283 


38 39 


(30,8) ^ (32,7) 


2.0283 


2.4396 


39 ^ 40 


(32,7) (34,6) 


0.2223 


0.4179 


39 — > 40 


(31,9) (32,8) 


3.6120 


> 5.0000 


40 ^ 41 


(32,8) -> (34, 7) 


0.4179 


0.5521 


41 42 


(32,9) (34,8) 


0.8329 


1.1372 


42 43 


(32,10)^(34,9) 


1.8473 


3.6391 


49 50 


(36,12,1) ^ (38,12) 


1.3753 


1.5634 


55 56 


(42, 12,1) -.(44, 12) 


0.4964 


0.5150 


59 ^ 60 


(47, 12,1) -.(48, 12) 


0.0000 


0.1024 



Table 3. Table of anomalies of the 2nd kind. Left column shows the change of the total 
particle number by one and column two the associated configuration change. The third and 
fourth column give the range of screening parameters where this transition occurs. 

and for k < h. Since k = 5 corresponds to a pair interaction of very short-range one might 
expect that further increase of k will not change the cluster structure qualitatively. To verify 
whether this is the case we performed, for all A^, additional calculations for an even larger 
screening, k = 20.0, cf. Fig. HI In most cases there is, indeed, no further change of the 
ground state configuration compared to k = 5, as expected. For three particle numbers, 
N = 44, 57, 58, the ground state configuration still changes in the "normal" way such that 
one particle is relocated from the outer shell to the inner shell. 

However, there are six remarkable cases which violate this trend: = 35, 36, 37, 39, 54 
and = 43. Consider first the total particle number = 43. This case is interesting because 
it is the only case where the central configuration (11, 1) is part of the ground state, apart 
from the cluster with 12 at screenings k > 4.1. This arrangement does otherwise not occur 
because the clusters prefer the platonic body with 12 particles on the inner shell (closed shell 
configuration). Here the configuration (31, 11, 1) becomes the ground state at k = 17.4 and 
remains the ground state for larger screening. 

In the other five cases we observed, at k = 20.0, several stationary states which differed 
only very little in their energies. We, therefore, performed molecular dynamics simulations 
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conf. 1 




conf. 2 




conf. 3 


35 


(29,6) 


3.2 


(28, 7) 


6.84 


(29,6) 


36 


(30,6) 


3.6 


(29,7) 


8.76 


(30,6) 


37 


(30,7) 


4.2 


(29,8) 


6.91 


(30,7) 


39 


(31,8) 


3.6 


(30,9) 


13.40 


(31,8) 


54 


(41,12,1) 


1.68 


(40, 13, 1) 


5.04 








5.04 


(39, 14, 1) 


15.02 


(41,12,1) 



Table 4. Table of anomalies of the 3rd kind. The 1st configuration is the ground state 
configuration up to the critical screening n^ri, then the ground state configuration changes 
in the standard way by adding a particle on a inner shell. This configuration then is the 
ground state up to the critical screening Kcra^ at which the cluster changes its ground state 
configuration back to the one it had a lower screening {anomalie of the 3rd kind). 

using separately each of these states as an input at k = 20.0 and then decreased the screening 
slightly, letting the system relax into a new stationary state, often with the same configuration 
and symmetry. This way we could be certain to follow all metastable states and independently 
record their energy dependence on k. The above five cases fall into two groups which differ 
with respect to the cluster symmetry. For the first, i.e. = 35, 36, 37, 39, the cluster 
decreases the number of particles on the inner shell when the screening is increased between 
K = 5.0 to 20.0. The resulting new ground state configuration contains again a platonic body 
on the inner shell. Allowing for such a highly symmetric configuration here turns out to be 
energetically more favorable compared to the previous shell configuration or a simple increase 
of the number of particles on the inner shell. 

Consider, for example, the cluster N = 39, cf. Fig.[6ld. Here, at k = 5 the ground state 
is (30, 9) until, at Kcr = 13.40, the configuration (31, 8) with one particle less on the inner 
shell becomes the ground state. 

The particle number reduction on inner shells is sometimes accompanied by another trend: 
with increasing n shells tend to split into subshells with close radii, as was already observed 
in Refs. fllTl l23l . This is observed e.g. for N = 35, cf. Fig. [6la. Here the configuration 
(28, 7) which is the ground state at k = 5 has in fact two subshells each containing 14 
particles which we will denote ([14, 14], 7). The radii of the two subshells differ only 
slightly, R2,i = 0.947 and i?2,2 = 0.875, respectively, while the inner shell radius is 
Ri = 0.426, clearly distinguishable from the outer shell. At K^r = 6.84 we observe a 
transition ([14, 14], 7) — > ([18, 11], 6), i.e. one particle from the inner shell moves outward 
and, in addition, three particles from the inner subshell move to the outer subshell. 

Similar behavior is observed for = 36, cf. Fig. [6lb. Here the configuration ( [15, 13] , 8) 
is the ground state at k = 5. At Kcr = 6.84 we observe a transition ([15, 13], 8) — ^ 
([22, 8], 6), where the inner shell loses two particles and, in addition, the inner subshell 
transfers 5 particles to the outer subshell. Analogously, for A^ = 37, cf. Fig. |6lc, the 
configuration ([12, 2, 15], 8) is the ground state at k = 5. At = 6.91 we observe a transition 
([12, 2, 15], 8) — > ([14, 16], 7), where the inner shell loses one particle. 

The cluster with A^ = 54 shows a similar behavior, cf. Fig. |71 Here, first the second 
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(29,6) , 



LU ^ 

b 



LU p 



' 1 ' (■22,8;'6)-(l5,13;8) — 




b) 

Kcr, = 8.76 
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K„^ = 13.40 
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r' (30,9) 


(31,8) 



4 6 8 10 12 14 16 18 20 
K K 

Figure 6. (colour online) Reentrant shell configuration changes for = 35 (top left), = 36 
(top right), iV = 37 (bottom left) and = 39 (bottom right). When k is increased, at Kcr i one 
particle moves towards the center and, at Kcr2, one particle returns to the outer shell restoring 
the former ground state configuration. The critical values of k are indicated by the vertical 
dashed lines. The solid red line shows the energy difference of the two configurations which 
has two zeroes. The legend shows the shell configurations, including the splitting of the outer 
sheU into subshells, given by the numbers in square brackets. 



shell popluation increases by one, at Hcri = 5.04, according to (40, 13, 1) — > (39, 14, 1). 
Further increase of screening makes a third configuration more favorable which has even 
two particles less on the second shell: the configuration (41, 12, 1) becomes the ground state 
again at Kcr2 = 15.02 which again is a consequence of the high symmetry (closed shell 
configuration). 

Finally, particularly interesting behavior is observed for all mentioned N = 
35, 36, 37, 39, 54, if a larger range of screening is considered, cf. Fig. [6] and Fig. |71 Here, 
for K > 1.7, there are two states which become the ground state. At a first critical value 
Kcri one particle moves to the inner shell until at Kcr2 this transition is reversed (the only 
difference for = 54 is that there is an additional ground state configuration between these 
two critical screenings): one particle moves outward and the original configuration with fewer 
particles on the inner shell is restored which remains the ground state for all larger values 
of K. This contradiction to the general trend (of increasing the inner shell population with 
increased screening), together with the reappearance of a ground state configuration, will be 
called "anomaly of the third kind". The complete set of these cases can be found in Tab. |4] 
with the exact critical screening parameters for the ground state configuration changes. The 
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Figure 7. (colour online) The energy differences of all states which become the ground state 
in the screening range 4.0 < k < 20.0 for the particle number N = 54. In this case increase 
of K leads to a ground state with fewer particles on the second shell, returning to a ground 
statea configuration, that already existed at lower screenings. In this case an "Anomaly of the 
third kind" is observed. 



reentrance of these ground state configurations at large screening are in all cases not different 
in their symmetry compared to the ground state configurations below Kcri, they have the same 
number of nearest neighbors and same shape of the Voronoi cells, with only their length scale 
strongly reduced due to the weaker interaction force. 

In general, the restored ground state configurations consist of platonic bodies on the inner 
shell, except for the case of = 37. Here, the ground state configuration changes from (30, 7) 
to (29, 8) at a screening value of k = 4.2 and back to (30, 7) at a screening value of k = 6.91. 
The 7 particles on the inner shell are not a platonic body and, as one can see from Fig. |4l it 
is not a common configuration compared to 6 or 8 particles in the center. Nevertheless, this 
can be understood by looking at the outer shell. The 30 particles are placed on the edges of 
an icosaeder which results in a highly symmetric configuration for the outer shell. 

4. Conclusions 

The goal of this paper was to present, for the first time, a detailed analysis of the ground state 
shell configurations of Yukawa clusters in a parabolic shperical confinement over a broad 
range of screening parameters k < 20. This allowed us to analyze the structural transitions 
occuring when the pair interaction changes from long range, in the Coulomb case, to short 
range, at the largest values of n. Focusing on a finite range of particle numbers, 11 < < 60 
we presented a complete overview on all existing changes of the shell configurations for 
K < 5. For larger k we also noted the cases where the configurations are different at k = 5 
and K = 20 (we cannot rule out that, in this range, there occur transtions in addition to those 
given). The general trend observed earlier [22] was confirmed: with increased screening, 
more particles tend to populate the inner shell(s) of the cluster giving rise to an average density 
profile which decreases increasingly fast towards the edge [|24ll25]| . There are, however, three 
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non-trivial deviations ("anomalies") from this shell filling sequence which were analyzed: 1.) 
upon K increase two particles move to (one of) the inner shell(s) at once. 2.) when the particle 
number is increased by one, at a fixed k, one particle moves from the inner to the outer shell 
and 3.) at very large k there exist cases of reentrent shell fillings: one particle returns from the 
inner to the outer shell. These anomalies are, in most cases, dictated by symmetry properties 
of the corresponding state which allow to lower the total energy. 

Our results are expected to be useful also for experiments with dusty plasmas and 
allow us to predict interesting parameter ranges which give information on the effect of 
symmetry on the structure of mesoscopic systems. In current experiments on spherical dust 
crystals performed at the Universities Kiel and Greifswald lfT3ll typical values of the screening 
parameter are in the range of 0.6 < k < 1.6. While this gives access only to a small 
part of the analyzed parameters where no reentrant shell fillings (third anomaly) occur, still 
the first two effects should be observable. While the experiments on small clusters do not 
necessarily reveal the ground state configurations, since often metastable states occur with a 
higher probability ll32l[33l the prediction of parameters where two states have the same energy 
is of practical interest for the analysis of potential energy barriers and intershell transitions. 
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